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Abstract. In this paper, we address a simplified version of a problem arising from vol¬ 
canology. Specifically, as reduced form of the boundary value problem for the Lame 
system, we consider a Neumann problem for harmonic functions in the half-space with 
a cavity C. Zero normal derivative is assumed at the boundary of the half-space; differ¬ 
ently, at dC, the normal derivative of the function is required to be given by an external 
datum g, corresponding to a pressure term exerted on the medium at dC. Under the 
assumption that the (pressurized) cavity is small with respect to the distance from the 
boundary of the half-space, we establish an asymptotic formula for the solution of the 
problem. Main ingredients are integral equation formulations of the harmonic solution 
of the Neumann problem and a spectral analysis of the integral operators involved in the 
problem. In the special case of a datum g which describes a constant pressure at dC, 
we recover a simplified representation based on a polarization tensor. 

Keywords. Asymptotic expansions; harmonic functions in the half-space; single and 
double layer potentials. 
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1. Introduction 

The aim of the paper is to provide a detailed mathematical study of a simplified version 
of a problem arising from volcanology. The analysis can be considered as a blueprint, 
useful to address the original problem in a forthcoming research; at the same time, the 
result has an interest on its own, entering in the stream of the asymptotic analysis for 
the conductivity equation, see [1, 2, 3, 4, 10] and references therein, with the principal 
novelties of dealing with a case of an unbounded domain with unbounded boundary and 
of a different choice of boundary datum (homogeneous on the boundary of the half-space 
and non-homogeneous on the boundary of the cavity). 

The geological problem is the detection of geometrical and physical features of magmatic 
reservoirs from changes within calderas. The starting evidence is that the magma exerts 

^Dipartimento di Matematica “G. Castelnuovo”, Sapienza - Universita di Roma, P.le Aldo Moro, 2 - 
00185 Roma (ITALY), aspriSmat.uniromal.it 

^Dipartimento di Matematica, Politecnico di Milano, Via Edoardo Bonardi - 20133 Milano (ITALY), 

elena.berettaSpolimi.it 

^Dipartimento di Matematica “G. Gastelnuovo”, Sapienza - Universita di Roma, P.le Aldo Moro, 
2 - 00185 Roma (ITALY), masciaOmat.imiromai.it AND Istituto per le Applicazioni del Galcolo, Gonsiglio 
Nazionale delle Ricerche (associated in the framework of the program “Intracellular Signalling”) 

1 



2 


a force aside the surrounding crust when migrates toward the earth’s surface, producing 
appreciable horizontal and vertical ground displacements, which can be detected by a 
variety of modern techniques (see the detailed description provided in [17]). As stated in 
[ 6 ], the main questions that emerge when monitoring voleanoes are how to eonstrain the 
souree of unrest, that is to estimate the parameters related to its depth, dimension, volume 
and pressure, and how to better asses hazards assoeiated with the unrest. For this purpose, 
the measurements of crust deformations are a useful tool to study magmatic processes 
since they are sensitive to changes in the source pressure and volume. 

From a modelling point of view, the displacements are described using the theory of 
linear elasticity, by replacing caldera with a half-space having a stress-free flat boundary, 
and the magma chamber by a cavity subjected to an internal pressure; in more detail, let 
V be the displacement vector v{x) = (ui(x), U2(x), U3(x)) and Vu = (Vu -|- (Vu)^)/2 the 
strain tensor, then the elastic model is defined by the linear system of equations 

(1) div(CVu) =0 in \ C, 


( 2 ) 


with boundary conditions 

(CVu)n ■ n = g on dC (CVu)n • r = 0 on dC 
(CVu)n = 0 on v{x) —)• 0 as jxj —)■ -|-oo, 

where C is the pressurized cavity, g is a vector-valued function represents the pressure, n 
the unit outer normal vector and r the tangential one; C is the isotropic elasticity tensor 
with Lame parameters A and ^ defined as 


C A /3 <Si I3 + 2/iIsym, 

with I 3 the identity matrix of order 3 and Igym the fourth order tensor such that Isym^ = A. 
The goal is to determine g and geometrical features of C (position, volume...) from surface 
measurements of the displacement field v. 

In applications, to handle more easily the inverse problem related to (l)-(2), the function 
g is taken constant and equal to a vector p G Additionally, from a geological point of 
view, sometimes it is reasonable to consider the magma chamber small compared to the 
distance from the boundary of the half-space, see [ 6 , 15, 17]; by adding these hypotheses 
and fixing the geometry of the cavity, some efforts have been done during the last decades 
to find some explicit or approximate solutions to the mathematical model. The simplest 
approximate solution obtained by asymptotic expansions is due to McTigue when the 
cavity is a sphere, see [15]. The other few solutions available concern ellipsoidal shapes, 
dike and faults, see [ 6 , 17]. 
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With the ultimate goal to study in detail the elastic problem (l)-(2), establishing an 
asymptotic expansion in the presence of a pressurized cavity of generic shape, in this paper 
we analyse a simplified scalar version of this model so as to shed light and mark the path to 
treat the elastic case. Denoting by Ml the half-space and M'^“^ its boundary, we consider 
the Laplace equation 

Art = 0 in Ml \ Ce 

with boundary conditions 

= q on dCe, —— =0 on M'^“^, u(x) —)■ 0 as |x| —-|-oo 

dn ^ ' dxd > V ^ II 

where Cg is the analogous of the pressurized cavity in the elastic case, with e a small 
parameter controlling its size, 5 is a function defined on dC^ and d > 3. Obviously, the 
choice to focus the attention on dimensions greater than two comes from the application 
we have in mind. 

Presentation of the main result. The main goals of this paper are first to analyse 
the well-posedness of the scalar problem, find a representation formula of the solution and 
then determine an asymptotic expansion with respect to the parameter e of the solution. 
For the first two steps we do not need to assume the cavity to be small. 

It is worth noticing that in terms of well-posedness the case of the half-space and, 
in general, of unbounded domain with unbounded boundary, is more difficult to treat 
compared to bounded or exterior domains since both the control of the solution decay and 
integrability on the boundary are needed. Indeed, it is typical to treat these problems 
by means of weighted Sobolev spaces, see for example [5]. In our case, we choose to use 
the particular symmetry of the half-space to prove the well-posedness in order to mantain 
a simple mathematical interpretation of the results. Therefore, we bring the problem 
back to an exterior domain in the whole space for which there is a vast literature on the 
well-posedness, see [9]. 

To trackle the issue of the asymptotic expansion, we consider the approach developed 
by Ammari and Kang based on single and double layer potentials for harmonic functions, 
see for example [3, 4]. This is the reason why we search an integral representation formula 
of the solution. To do this, we take advantage of the explicit expression of the Neumann 
function for the half-space 

N{x, y) = r{x -y) + r{x - y), 

where T is the fundamental solution of the Laplacian and x is the symmetric point of 
X with respect to the plane, in order to get a representation formula containing only 
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integral contributions on the boundary of C. In detail, we find that 



N{x,y)g{y) 


JL 

dn. 


N{x,y)f{y) 


da{y), 


X e Ml \ C, 


where / is the trace of the solution u on dC. From the point of view of the inverse 
problem we are interested in evaluating the solution u on the boundary of the half-space; 
since r{x — y) = r{x — y), for x G the integral formula becomes 


1 

2 



r{x 


y)g{y) 


d 

dUy 


r{x 


y)f{y) 


da{y), 


x G 


Taking B a bounded Lipschitz domain containing the origin and z G Ml we consider 
:= C = z + eB with the assumption that dist(z, M'^“^) > do > 0. Therefore, defining 
e) = g{z -|- eC), with C ^ B, and Ssg the single layer potential, the main result holds 


Theorem 1.1. For any e > 0, let g ^ L?‘{dCe) such that g is independent on e. Then, at 
any x G we have 

Ue{x) = 2e'^“^r(x - z) [ g{C)dcr{C) 

JdB 

+ 2e''vr(x -z)- {n^ 5 b5(C) - C5(C)} 

where denotes a quantity uniformly bounded by with C = C'(do) which tends 

to infinity when do goes to zero. 


Finally, with the asymptotic expansion in hand, we consider the Neumann boundary 
datum g = —p ■ n where p is a constant vector in M'^. This particular choice has a double 
purpose: to reconnect this problem with the constant boundary conditions of the elastic 
model and to make more explicit the integrals in the asymptotic formula. The result we 
get contains the same polarization tensor obtained by Friedman and Vogelius in [10] for 
cavities in a bounded domain. 

The organization of the paper is the following. Section 2 is divided into three parts: 
in the first one we recall some well-known results about harmonic functions and layer 
potentials; in the second one we examine the well-posedness of the scalar problem; in the 
third one we get the representation formula of the solution. In Section 3, we state and 
prove a spectral result crucial for the derivation of our main asymptotic expansion. In 
Section 4 we present and prove the theorem on the asymptotic expansion and finally we 
illustrate the result for the particular choice g = —p ■ n. 
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Notation. All the analysis is performed in with d > 3; Br{x) C denotes the ball with 
centre x and radius r > 0, and ujd the area of the (d — l)-dimensional unit sphere. The scalar 
product between two vectors is represented hy x ■ y and rix indicates the unit outward normal 
vector in x on the boundary of some specified domain. The fundamental solution T of the Laplace 
operator in with d > 3, is given by r(a::) = Kd\x\‘^~‘^ with := l/u}d{2 — d). 

Points X S d > 3, are decomposed as a; = {x',Xd) where x' = (a;i,- • • ,Xd-i)- We denote 
by Ml the half-space {x S M'^ : Xd < 0} and by M"^”^ its boundary. Given a point x G Ml, its 
reflected point {x', —Xd) with respect to the plane Xd = 0 is represented by x. 


2. The direct problem 


In this Section, we analyse the boundary value problem 


( 3 ) 




Au = 0 
du 
dn 
du 
dxd 
u —)• 0 


= 9 
= 0 


in Ml \ C 
on dC 
on M'^“^ 
as \x\ —)■ -boo 


where C is the cavity, with a twofold aim: to establish well-posedness of the problem and 
to provide a representation formula. 


Preliminaries. The specific symmetry of the half-space permits to show well-posedness 
by extending the problem to an exterior domain, viz. the complementary set of a bounded 
set. Hence, it is useful to recall the classical results on the asymptotic behaviour of 
harmonic functions in exterior domains. Given a bounded domain 1? C M*^ with d > 3, if u 
is harmonic in M'^ \ 17 then v is harmonic at infinity if and only if v tends to 0 as |x| —)■ oo. 
Moreover, there exist rg, M > 0, such that if |x| > rg, the following estimates hold 

(4) \v{x)\ < , \vx^{x)\ < M\x\^~‘^ <M\x\~'^ 


for i, j. A; = 1,..., d. The proof can be found in [9] (see also [8, 14]). 

The representation formula makes use of layer potentials whose definitions we now recall; 
see [3, 9, 14]. Given a bounded Lipschitz domain 17 C M'^ and a function (p G L^(5l7), we 
introduce the integral operators 

Sa(p{x) := [ r{x-y)if{y)da{y), x G M"*, 

Jon 

f 9 

Daip{x):= —r{x - y)ip{y)da{y), 

Jan duy 


( 5 ) 


X G M'^ \ dQ, 
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which are called, respectively, single and double layer potential relative to the set fi. By 
definition, Sqlp and Dfpp are well-defined and harmonic in \ dfi. Further, still for 
d> 3, we have 

= 0 {\x\‘^~'^) and D^ip = 0 {\x\^~'^) 

as |x| — 7 - -|-oo. In addition, if p has zero mean on df2, the decay rate of the single layer 
potential Sq is increased, precisely, 

if [ ip{x)da{x) = 0, then Sqp = 0 {\x\^~’^) 

JdQ 

as \x\ —)■ -|-oo (such property holds for d > 2). 

Next, we introduce the compact operator Kq : L^^dQ) —?■ L^(cB7) 

( 6 ) Kq(p{x) := —p.v. f — — ^iy)dcr{y) 

i^d Jan \x - 2/r 

and its L^—adjoint 

(7) K} 2 ^{x) = —p.v. f ^ — '^^^—^ip{y)da{y). 

Wd JdQ \x-y\^ 

Given a function v defined in a neighbourhood of <917, set 


•= li™ v{x±hnx), X & df2. 


h^ 0 + 


The following relations hold, a.e. in df2 

( 8 ) SqpI = Sq(p\_, 


dux 


± 


[±ll + K},)p, D^p\ ={^^ll + Kn) 


if. 


For the proof see [3, 9] . 


Well-posedness. Proving existence and uniqueness results for unbounded domains with 
unbounded boundary is, in general, much more difficult with respect to the case of exterior 
domains. The main obstacle is the control of both solution decay and integrability on the 
boundary and usual approach is based on the use of weighted Sobolev spaces [5]. Here, 
we take advantage of the symmetry property of the half-space to extend the problem to 
the whole space and to establish well-posedness resorting in a standard Sobolev setting. 
Given a bounded Lipschitz domain C C Ml and the function g : dC —)• M, we define 

C := {{x', Xd) : (x', -Xd) G C} 

and G : dC U dC —)■ M as 


G(x): 


g{x) if X G dC 
g{x) if X G dC. 
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Theorem 2.1. The problem (3) has a unique solution. This solution coincides with the 
restriction to the half-space Ml of the solution to 

T 

(9) AU = 0 in M'^ \ U = G on dC U dC, U ^ 0 as |x| —)■ +oo. 


Proof. The proof is divided into three steps: uniqueness for (9), existence for (9), equiva¬ 
lence between (3) and (9). 

1. For T := CU C, let -R > 0 be such that A C and set f2ji := Bn \ A. Given two 

solutions, Ui and C/ 2 , to problem (9), the difference W := Ui — U2 solves the corresponding 
homogeneous problem. Multiplying equation AW = 0 by IF and integrating over the 
domain Qn = Bn \ A, we infer 


0= / W{x)AW{x)dx 
d 




W{x)^W{x)da{x) - \VW{x)\ dx, 

JdBnio) on Jaji 

using integration by parts and boundary conditions. Exploiting the behaviour of the 

harmonic functions in exterior domains, as described in (4), we get 


/ W{x)-^W{x) da{x) 

IdBnio) on 


< 


C 


Rd-2- 


Then, as i? —)• 00 , we find 


\VW{x)fdx = 0 


which implies IF = 0. 

2. We represent the solution of (9) by means of single layer potential 


( 10 ) 


SAfjix) = [ r{x-y)'if{y)da{y), x G M'^ \ T, 
JdA 


with function ip to be determined. By the properties of single layer potential, SAtf is 
harmonic in M*^ \ dA, S a'4>{x) = 0{\x\‘^~^) as |x| —00 and we have 


dSAfj 

dn 


(x) 


= \'p + 


We now prove that there exists a function if such that 


X G dA. 


(11) (^I + 'p{x) = G{x)., X G (9/1. 

To this aim we prove that the operator ^I+K\ is injective. Indeed, given Q G 
define V := SaC- Then, from the properties of layer potentials, VF solves 

d 

AF = 0 in M*^ \ T, —F = 0 on 9/1, F —)• 0 as |x| —)■ 00 , 

on 







8 


hence, from Step 1., the function V is identically zero; then it follows that C = 0. Finally, 
observing that K\ is a compact operator and G G L‘^{dA), equation (11) admits a unique 
solution. 

3. To show that u := we have to verify that the normal derivative is null on the 

boundary of the half-space. This is an immediate consequence of the symmetry property 

(12) U{x',-Xd) = U{x',Xd), 


which follows from the observation that, by definition of the boundary datum G, the 
function u{x',Xd) ■= U{x', —Xd) solves (9) and the solution to such problem is unique. 

As a consequence of (12), we obtain 

du , dU , dU , 

Thus the derivative of U with respect to Xd computed at any point with = 0 is zero. □ 


Representation formula. Next, we derive an integral representation formula for the 
solution u to problem (3). This makes use of the single and double layer potentials defined 
in (5) and of contributions relative to the image cavity C, given by 

Sc‘p{x) ■■= f r{x -y)(f{y)da{y), xGR'^, 

~ T d 

Dc(p{x) := / -7^r{x - y)(p{y)da{y) x G \ dG. 

JdC ony 

These operators, referred to as image layer potentials, can be read as single and dou¬ 
ble layer potentials on G applied to the reflection of the function ip with respect to Xd 
coordinate. 


Theorem 2.2. The solution u to problem (3) is sueh that 

(14) u{x) = Scg{x) - Dcf{x)+ Scg{x) - Dcf{x), x G Ml \ C, 

where Sc,Dc are defined in (5), Sc,Dc in (13), g is the Neumann boundary eondition 
in (3) and f is the traee of u on dG. 

Using properties of layer potentials, from equation (14), we infer 

/(x) = Scg{x) - [-\l + Kc^ f{x) - Dcf{x) + Scg{x), x G dG, 

where Kc is defined in (6). Thus, the trace / satisfies the relation 

( 5 / -h Kc + Dc^ f = Scg + Scg, 

which will turn out to be useful in the sequel. 
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Before proving Theorem 2.2, we first recall the definition of the Neumann function, see 
[11], that is the solution N = N(x,y) to 

AyN{x,y) = 6x{y) in Ml, ^(x,?/) = 0 on M'^“\ 

where 6 x{y) is the delta function centred in a fixed point x £ M'^ and dN/dyd represents 
the normal derivative on the boundary of the half-space Ml. The classical method of 
images provides the explicit expression 


N{x,y) 


Kd 


\x - y\ 


d-2 




l^d 


\x - y\ 


d-2 ■ 


With the function N at hand, the representation formula (14) can be equivalently written 
as 


u{x) =J\f{f,g){x) 


(15) 


IdC 


d 


N{x,y)g{y) - —N{x,y)f{y) 

OUv 


da{y), X £ Ml \ C, 


which we now prove. 


Proof of Theorem 2.2. Given i?, e > 0 such that C C Bij(O) and B^{x) C Ml \ C, let 

■= (m 1 n Br( 0)) \ (c U B,{x)). 

We also define dB^{0) as the intersection of the hemisphere with the boundary of the 
half-space, and with dB’^{0) the spherical cap (see Figure 1). Applying second Green’s 



Figure 1. Domain used to obtain the integral representation formula (14). 
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identity to N{x, •) and u in we get 


0 = 


[N{x, y)Au{y) - u{y)AyN{x, y)] dy 


+ 


^R, 

f 

dB^iO) 

f 


f)qj f) 

N{x,y) — {y) - —N{x,y)u{y) 


f)qj f) 

N{x,y) — {y)-—Nix,y)u{y) 


dm. 


dm, 


+ 


/ dB£{x) 


d dvj 

N{x,y)u{y) - N{x,y)^—{y) 


dn. 


dm, 


da{y) 
da{y) 
da{y) 


lac 




da{y) 


=Ii +l 2 + h -M{f,g){x). 

The term Ii is zero since both the normal derivative of the function N and u are zero 
above the boundary of the half-space. 

Next, taking into account the behaviour of harmonic functions in exterior domains, 
formulas (4), we deduce 


d 


IdBUo) duy 


N{x,y)u{y) da{y) 


IdBUo) 


N{x,y)^^^^da{y) 


< 


< 


C 

C 

i?2d-3 


/ , da{y) = 

/ ^ da{y) = , 

JdBt(O) R 


C 

c 


where C denotes a generic positive constant. As i? —>■ -|-oo, I 2 tends to zero. 
Finally, we decompose I 3 as 


h — hi — I 32 — 


d 


N{x,y)u{y) da{y) - 


ldBs{x) dfly 

Using the expression of N and the continuity of u, we derive 


f N{x,y)^{y)da{y). 

dB,{x) dUy 


dsi = 


IdBeix) OUy 


d_ 

dm 


N{x,y)u{y)da{y) =u(x) 


f - 

ldBs{x) dlfly 


+ 


N{x,y)da{y) 
d 


' dBs{x) 


u{y) -u{x)] ^—N{x,y)da{y), 
dn. 


which tends to u{x) as e —)• 0. Moreover, we infer 

du 


\h 2 \ < C sup 

y£dBe{x) 

< C' sup 

y£dBs(x) 


dUy 

du 




dm. 


ldB,(x) £' 


N{x,y) da{y) 

^da{y) -h [ 

‘ Ja 


aB,ix) \x - y\^ 2 


CJ?da{y) 
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Observing that both the integrals tend to zero when e goes to zero because the second 
one has a continuous kernel while the first one behaves as 0 (e), we infer that 1^2 —t 0 as 
e — 7 - 0. Putting together all the results, we obtain (15). □ 

3. Spectral analysis 

Following the approach of Ammari and Kang, see [3, 4], in this Section, we prove the 
invertibility of the operator \ l + Kc + Dc showing that, under suitable assumptions, the 
following inclusion holds 

a{Kc + Dc) C (-1/2,1/2], 

Such task is accomplished by determining the spectrum of the adjoint operator Kq + 
in L?‘{dC), relying on the fact that the two spectra are conjugate. 

The explicit expression of is in (7). Computing the L^-adjoint of Dc is straightfor¬ 
ward: indeed, given G L?‘{dC), we have 

[ 'il>{x)Dcip{x)da{x) = ( ijjix) (— f da{x) 

Jac Jac \aJdJac \x-yr ) 

= / < 7 >(y) (— / il;{x)da{x)] da{y) 

Jac V<^d Jac k - yp J 

and thus 

(16) D*c<f{x) = — f (^ '^^ ^{y)da{y). 

dd Jac \y - x\^ 

Note that the kernel of the integral operator is smooth on dC. 

As proved in [13], for smooth domains the eigenvalues of on L?‘{dC) lie in (—1/2,1/2]; 
the same result it is also true for Lipschitz domain (see [3, 7]). With the same approach, 
it can be shown that the same property holds true for + D^. 

Theorem 3.1. Let C he an open bounded domain with Lipsehitz boundary. Then 

a{K*c + D*c)ci-l/2,l/2]. 

For completeness, we provide here a complete proof of such fact. 

Firstly, we observe that the regular operator Dq on the boundary of the cavity can be 
seen as the normal derivative of an appropriate single layer potential. 

Lemma 3.2. Given ip G L^{dC) we have that 

where ip G L‘^{dC) is defined by ip{x) := ip{x). 


x G dC, 
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Proof. Using the expression (16) of Dq and the identity 


V. 


{2 — d)\x — J |x —y|‘^’ 


x-y 


we find that 


D*(jif{x) = Vj 


'ac \y - X 


i^d p{y) , r \ \ 
—^da{y) ) • 


Given ip G L?‘{dC) and ip G L‘^{dC) as previously defined, we have 


Jdc\y-x\'^^ Jac \z - xl^^-"^ 


ip{z) 


lac \y - x\^- 


Lc = L 


ip{z) 


which gives the conclusion. 

We are now ready to prove the main result of this Section. 


□ 


Proof of Theorem (3.1). Given ip G L'^{dC), let if be defined by if '■= ScT + Sp;ip. By the 
known properties of single layer potentials, we derive on dC 


dip 

dn 


= (±\I + Kfj + D}j)ip 


and, as a consequence. 


(17) 


dip 

dn 


+ 


dip 

dn 


= 2(k*c + D*c) ip, 


dip 

dn 


dip 

dn 


= T- 


Taking a linear combination of the two relations in (17), we deduce 


\ X f 


dip 

\ 1 f dip 

dpJ 


+ 

dn 

2 ydn 

dn 


= (A- 1 ^ ^ 
' 2 / dn 




If A is an eigenvalue of Kq + Df; with eigenfunction ip, then 


A- 


1 \ dip 


2 J dn 


- A + 


1 \ dip 


2 J dn 


= 0 , 


on dC. 


Multiplying such relation by the function ip and integrating over dC, we get 
1 


(18) 


A- 


lac dn 


da{x)-{x+\) 


da{x) = 0. 
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Integrating by parts we have 


(19) 


/ daix) = [ 'ip{x)Ail;(x)dx+ [ IVipix)]^ dx 

Jac on _ Jc Jc 

= [ |Vi/i(x)Pdx. 

Jc 


The first integral in (18) can be dealt with as done in the proof of Theorem 2.2. Precisely, 
given large i? > 0, applying the Green’s formula in fJn := n 5^(0)^ \ C, we get 


[ 'ip{x)^{x) da{x) = [ ip{x)^{x)da{x) + [ ip{x)^{x) 

Jac on , JdB>X(0) OXd JdBUo) On 


dxp 


d'ljj 


da{x) 




'0(x)A-0(x)dx — / |VV'(x)| dx, 


' Qt. 


where is the intersection of the hemisphere with the half-space and dB^{0) is the 

spherical cap. The quantity d'ljj jdn is identically zero on the boundary of the half-space 
since the kernel of the operator is the normal derivative of the Neumann function which, 
by hypothesis, is null on Moreover, ij: is harmonic in Qr, so we infer 

d'ljj 

ipyx] 

Recalling the asymptotic behaviour of simple layer potential, 

l'S'c‘^1 + = 0{\x\^~'^), IV^cV^I + |V5g;(/?| = 0(|x|^“‘’*) as \x\ —oo. 

we obtain, for some C > 0, 


/ ^{x)^{x) 
lac dn 


da{x) = / i/)(x)^(x) da{x) — [ \V'tp{x)f dx. 

JdBi(o) dn , Jn„ 


IdBUo) dn 


da{x] 


< 


< 


JdBj^iO) 

C 

R2rf-3 ^ 




da{x) 

f da{x) = 
obUo) R 


Passing to the limit R —)■ -|-oo, we find 

( 20 ) ‘ 


[ 'ijj{x)— da{x) = — [ _\V'ijj{x)\‘^dx. 

Jac dn Jr^xc 

Plugging (19) and (20) into (18), we infer the identity 

A - / |Vi/)(x)Pdx-I-(^A-I-/ |VV’(x)Pdx = 0, 

2/ 7r^\c \ 2 / Jc 


that is 

with 


{A + B)X=^{A-B) 


A := 


:= / \V'ijj{x)\^dx and R := / |VV’(x)l^dx. 

jR-jXC Jc 
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The coefficient of A is non-zero. On the contrary, \i A + B = {) then = 0 in Ml which 
means that = 0, hence, from the second equation in (17), we get </? = 0 in dC. 
Therefore, solving with respect to A, we finally get 


( 21 ) 


A 


1 A-B ^ 

2 ‘ A + B ^ 


■ 1 1 ' 
~ 2 ' 2 


The value A = —1/2 is not an eigenvalue for the operator Kq + D^. Indeed, in such a 
case, we would have 

A = [ |Vi/(x)P dx = 0, 

Jr'^\c 

and thus ?/ = 0 in Ml \ C. By definition of if:, we deduce that 'i/ = 0 on dC and since ijj 
is harmonic in C, we get that = 0 also in C. As before, by (17), this would imply that 
(/? = 0 in dC. □ 


For completeness, let us observe that the value A = 1/2 is an eigenvalue with geometric 
multiplicity equal to one. Indeed, identity (21) implies that, for such value of A, 

B = [ dx = 0, 

Jc 

hence "i/ is constant in C. Normalizing ip = 1 in C, the function ip in Ml \ (7 is given by the 
restriction of the solution U to the Dirichlet problem in the exterior domain M'^ \ (^C U 
with boundary data equal to 1. Then, by the second equation in (17), the function ip is 
the normal derivative of U at dC. 


4. Asymptotic expansion 

In this Section, we derive an asymptotic formula for the solution of problem (3) when 
the cavity C is small compared to the distance from the half-space For the reader’s 

convenience, we recall that the cavity C has the structure 

Ce := C = z + sB 

where B is a bounded Lipschitz set containing the origin. Moreover, we assume that 

(22) dist(2;,M'^“^) > do > 0 

otherwise, for the application we have in mind, the problem does not have a real physical 
meaning. To emphasize the dependence of the solution to the direct problem by the 
parameter s we denote it by Ue- For brevity, we denote the layer potentials relative to Cg 
by the index s, viz. 

S, = Sc,, De = Dc,, Se = Sc,, De = Dc,, Ke = Kc, 





15 


and the trace of the solution on dC^ by f^- In this way the representation formula (14) 
reads as 

Us = SeQ - Defe “ Defe + SeQ- 

At X G taking into account that x = x, it follows that 

Seg{x)= [ r{x-y)g{y)do{y)= j r{x - y)g{y) da{y) = Ssg{x) 

JdCe JdC^ 

Dsfeix) = f ^r{x - y)fsiy) da{y) = [ -^^{x - y)fs{y) da{y) = Defs{x) 

J dCg J dCs: C/'^y 

Hence, we obtain the equality 

5 Ue{x) = Seg{x) - DJs{x), X G 
Associating with the relation at the boundary dCs and by (3.1)] 

fe{x) = + Ke + De^ (^Ssg + Ssg'^ {x), x £ dCe, 

we get the identity 

(23) ^Usix) = Ji(x) + J2{x), X G 
where 

Ji(x) := / r{x - y)g{y)da{y), 

JdC^ 

J 2 {x) := - (5^ + Ks + (^S^g + Seg^ (y) da{y). 

Analyzing in details the dependence with respect to e of such relation, we obtain an explicit 
expression for the first two terms in the asymptotic expansion of at 

In what follows, for any fixed value of e > 0, given h : dC^ —)■ M, we introduce the 
function h : dB —)■ M defined by 

^(C; e) := h(2: + eC), (^dB. 

This definition is useful to consider integrals over a set that is independent on e. 

Theorem 4.1. Let us assume (22). For any e > 0, let g £ Li^^dC^) sueh that g is 
independent on e. Then, at any x £ the following expansion holds 

Ue{x) = 2 e‘^~^r{x - z) [ g{C)dcr{C) 

(24) 

+ 2e^Vr(x - [nc SegiC) - Cdio] da{C) + 0{e^+^), 

where denotes a quantity uniformly bounded by with C = C'(5o) whieh tends 

to infinity when do goes to zero. 





To prove the theorem we first show the following expansion for the operator 
Lemma 4.2. We have 

(25) (p + K, + (Ssg + S^g) (z + eC) = £ (p + Kb) SbMC) + 

Proof. We analyse, separetely, the terms ^p + + D^) and 8^ + Ss, collecting, at the 

very end, the corresponding expansions. 

Since K^ + De is compact and its spectrum is contained in (—1/2,1/2] there exists (5 > 0 
such that 

a{K, + D,) C (-1/2+ 5,1/2], 

Then, the operator 

A, := ll-K,- D, 

is such that c {As) C [0,1 — S) and thus has spectral radius strictly smaller than 1. As a 
consequence, taking the powers of the operator As one finds 


(26) jjA^II < 1 V/i and ||A(i°|| < 1 for some/ iq- 

The inverse operator of / — = p + Ks + Ds can be represented by the Neumann series 

+ 00 +00 

'll-Ks-D, 


that is 


{I - As)-^ = Y.A^, = - K, 


/i=0 


/i=0 


At the point z + ef, we obtain 


(^Ks + Ds) (p{z + eC) 


= —p.v 




Wd ■ ' JdCe \z + £C-y\ 


1 

= —p.v. 


dB 

= KBm + e^~'ReP{0, 


iv C) _g gci 1 


/ — 

'dc, duy 

• d 


r{z + eC - y)ip{y) da{y) 


IdB dUrj 


r{z + ef — z — £g)(p{g) da{g) 


where 


/* 

RePiC) ■= -^R (z-z + £{l-g)) (p{g)da{g) 


IdB duy 

is uniformly bounded in s. Using these results, we calculate A^ highlighting the term that 

do not contain e and the one of order d — 1, that is 

, h 

^a— I T . 


Ai = 






where 
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For ho as in (26) and h > ho we have 
where [ • ] denotes the integer part, and thus 

+00 +O0 

h=0 h=0 

giving the absolute convergence of J2Eh,e- Summarizing we conclude that 

(27) {I-Ae)-^ = [^I + KBy\o{e^-^). 

Let us evaluate the term Ss + Ss- We have 

Seg{z + eC)=[ r{z + eC - y)giy)da{y) 

JdC. 


= e 


IdB 


r{C-0)g{O)daie) = eSBg{C) 


and 


Seg{z + eC)= ri^z + eC-y) g{y)da{y) 

dCc 


= [ r ( 

JdB ^ 


z — z + e(C — 9)) g{d)do{9) 


= ^r{z — z) [ g{6)da{9) + 0{e^) 
JdB 

where we have used the zero order expansion for F. 

Collecting we infer 


[S,g + S,g) {z + eC) = eSBg{C) + 0{e^-^) 

and combining with (27) we obtain the conclusion. □ 

Proof of Theorem f.l. To prove (24), we analyse the two integrals Ji and J 2 in (23). 

For X, C G with x 7 ^ 0 and e sufficiently small, we have 

r(x - eC) = E{x) - e Vr(x) • C + 0{e^). 

Hence, for x G we get 

Ji = e'^~^f r{x - z - eC)g{C)dcr{C) 

(28) 

= e^^r{x-z) g{C)da{C) - e^Vr{x - z) ■ Cg{0d<r{0 + O{e'^^^). 

Job Job 

Next we consider the second integral in (23), written as 

J 2 = [ 7^F(x - z - eC)heiC) dcr{C), 

JdB OUq 
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where the function is given by 

(29) = ( 5 -^ + Ki; + {SeQ + SeQ^ {z + 6 ^) 

For X, C G with x 7 ^ 0 and e sufficiently small, there holds 

(30) Vxr{x + eC) =Vxr{x) + 0{e), 
therefore, taking advantage of the expansion (25), 

J 2 = e'^“^ [ -^r{x - z)he{C} da{C) + 0{£‘^) 

JdB onc^ 

= -^r{x - z) [\I + Kb) SugiC) da{0 + 0(e''+^). 

Collecting the expansions for Ji and J 2 , we deduce (24). □ 

We show that the term + Kb) SBg{x), for x G dB, represents the trace of the 
solution of the external domain related to the set B and with Neumann boundary condition 
given by g. To this aim, we consider the problem 

dU 

(31) AU = 0 mW^\B, = g on dB, U —>■ 0 as |x| —)■ + 00 , 

on 

where the cavity B is such that 0 € B. 

Proposition 4.3. Let us define /(x) := then 

+ SBg{x) = f{x). 

Proof. As done in the proof of Theorem 2.2, that is, applying second Green’s identity to 
the fundamental solution P and U in the domain Bb,{0) \ B, with R sufficiently large, it 
can be proven that the representation formula for U is 

(32) U{x) = SBg{x) - DBfix ), xGR<^\B. 

Therefore, from single and double layer potentials properties 

f{x) = SBg{x) - (^-^I + Kb'^ f{x), xedB, 

hence ^ 

f{x) = Ql + Kb^ SBg{x), X G dB, 

that is the assertion. □ 

Now, we want to consider a specific case of the Neumann condition on the boundary of 
the cavity so to get an explicit expression of the asymptotic expansion in terms of the 
polarization tensor and the fundamental solution. 
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Corollary 4.4. Given p G let the boundary datum given by 


g = —p ■ n. 

Then, the following expansion holds 

(33) Ue{x) = 2s'^Vr{x - z) ■ Mp + 0{e^+^), x G 

where M is the symmetrie positive definite tensor given by 


(34) 


M := / <8) (C +'I'(C)) f^cr(C) 

JdB 


and the auxiliary funetion I? has eomponents i = 1,... ,d, solving 


Aifi = 0 in M“ \ B, 
Proof. Let us set 


dn 


= —n,- on dB, 




0 as X —)■ + 00 . 


Ji := Vr(x - H Sb[-P ■ n](C) da{C), 

J2 := Vr(x - z) • / Cp-n(-da{C). 

JdB 

Then, expansion (24) with g = —p ■ n gives 

hueix) = -e'^“^r(x - z) [ p-nQ da{Cj + Ji + J 2 + 

(35) JdB 

= Ji + J 2 + 0{e^+^) 

since divergence Theorem guarantees that the first term in the expansion for is null. 

From the equation (31), with g = —p ■ n, since the problem for U is linear, we can 
decompose U as U = J2i Ui where the functions C/j, for i = 1, • • • ,d, solve 


dU- 

AUi = 0 in \ i3, = —piUi on dB, [/, 

on 


0 as |x| —+ 00 . 


From the definition of the functions ihi, we deduce U = p ■ . Using Proposition 4.3, the 

term Ji can be rewritten as 

Ji = Vr{x-z)-[ ('L(C) • p) do-(C) = Vr(x - z) • / (n^ (g) 4'(C))pdT(C). 

JdB JdB 

To deal with the term J 2 , we preliminarly observe that 

/ da{C) = f (nc;0C) da(C). 

JdB JdB 
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Indeed, for ... n(^^d), for any i,j G {1,..., d}, it follows 

/ Cin(^,jda{C)= nQ-Ciejda{C)= div {Ciej) dC = ej-CidC 
JOB JdB JB Jb 

= / V(Cj)-eidC= / div{Cjei)dC= / nc_ ■ Cjeida{C) 

JB JB JOB 


IdB 


ncj Cj da{C) 


where ej is the j-th unit vector of Thus, we get 

J2 = Vr{x - z) ■ [ {C'^nc)pda{C) = Vr{x - z) ■ [ {nc <Si C) pda{C). 

JdB JdB 

Collecting the expressions for Ji and J 2 , we obtain formula (33). 

Symmetry of the tensor M, defined in (34), follows from 

54' o 


/ ^i{C)nc,j daiC) = - [ d>i{C)-^{C) daiC) 
JdB JdB on 


div (4/^(0V4/,(C)) dC = 


V4/i(C)-V4/,(C)dC 


where the last term is obviously symmetric. Taking p G M'^, we consider 

p- Mp= [ {ne-p){C- p) da{C) + ( (n^ • r/)(4'(C) • p) da{C) = Ii + h- 
JdB JdB 

The positivity of the tensor follows from the divergence theorem and integration by parts, 
in fact 

h= f div{{C ■ p)p) dC = [ p-V{C-p)dC= f \p\'^ dC = \p\‘^\B\ 

Jb Jb Jb 

In the same way, by the definition of the function 4^ 


h = - 


f ■v)da{p) = f div((4'(C) •r/)V(4'(C)-r?)) do-(C) 

JdB on jRd\B 

= f |V(4/(C)-ry)|'dC 

jRd.\B 


The sum of Ii and I 2 gives the positivity. 


□ 


For specific forms of the cavity C, the auxiliary function 4 / can be determined explicitly, 
providing a corresponding explicit formula for the polarization tensor M. The basic case 
is the one of a spherical cavity (see [10]). If i? = {x G : jxj < 1}, then a direct 
calculation shows that, for i = 1,2,3, there holds 4'i(x) = Xi/{{d — l)|x|"'), and thus 

MC) = j^Ci, C&dB. 

As a consequence, the polarization tensor is a multiple of the identity and, precisely, 

M = -\B\I = 2ttI. 

2 ' ' 
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Then, the asymptotic expansion (33) becomes 

Ue{x) = 47re‘^Vr(x — z) ■ p + x G 

Explicit formulas can be provided also in the case of ellipsoidal cavities (see [1, 3, 4]). 

In general, for given shapes of the cavity C, such auxiliary function can be numerically 

approximated and, thus, the first term in the expansion (33) can be considered as known 

in practical cases. 
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